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Introduction 



In 1963 Staruszkiewicz considered for the first time gravity in 2+1 dimensions cou- 
pled to point-particles [Q. In 1984 the subject was reconsidered by Gott, Alpert, 
Giddings, Abbot, Kuchar, Deser, Jackiw and 't Hooft 0. In these articles they 
found solutions for the gravitational field around N static point-particles. They 
also solved the case of one spinning particle located at the origin. The key result 
was that locally these space-times are flat except at the particles positions. This 
implies that the gravitational field contains no degrees of freedom as can also be 
deduced from a simple counting argument: We have 6 independent metric com- 
ponents minus 3 first-class constraints minus 3 gauge-fixing conditions, resulting 
in degrees of freedom. The reason that these spaces are not completely trivial 
is because of non-trivial boundary conditions on the (flat) coordinates. For in- 
stance, in the case of a massive point-particle sitting at the rest in the origin, we 
have to cut a wedge out of space-time and identify opposite points of the wedge. 
So the range of the angular coordinate ip is: G [0, 27r(l — AmG)] (m is the mass 
of the particle). In general this identification condition on the coordinates is not 
a simple rotation but a Poincare-transformation (see section 1). Not much later 
Deser and Jackiw found solutions for gravity with a cosmological constant. The 
case A > (de Sitter space) coupled to 2 static, antipodal particles was solved 
[0. In a geometrical approach 't Hooft solved the N-particle case (A = 0), and 
proved that a Cauchy-formulation was possible within which no closed timelike 
curves could occur A different view on the problem was provided by Achu- 
carro and Townsend and later by Witten They considered a Chern- Simons 
theory with a gauge-field taking values in the Poincare-algebra, and proved 
that this theory is equivalent to 2+1 dimensional gravity. Later Grignani and 
Nardelli invented a consistent way to couple point-particles to this gauge field |^. 
This Chern-Simons approach is closely related to the description of gravity using 
Ashtekar variables It is however not known to me if people considered the 
coupling of point-particles (For a review on loop-quantisation one should consult 
the lecture by Kirill Krasnov in this volume). Finally I would like to mention 
Waelbroeck's approach [jl4[ who consideres these Ashtekar variables on a lattice 
in order to obtain a finite set of degrees of freedom. This is an exact description 
because only the handles and particles are the true degrees of freedom. 

A first step towards quantisation was made by Mazur, 't Hooft, Deser, Jackiw 
and de Sousa Gerbert who studied the scattering of 2 quantum particles. ||^. 
Later the calculation was also done in the Chern-Simons approach Carlip 
considered the scattering of N particles where he stressed the role of the braid- 
group in 2+1-gravity |jlO|. All these approaches have their own way of quantising 
the theory. The problem is however that not all these quantisations seem to be 
equivalent as shown by Carlip p . In a way this is disappointing, but it reflects the 
fact that the problems encountered in quantising 3+1-dimensional gravity still 
survive the dimensional reduction to 2+1 dimensions. As the theory contains 
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no gravitons, and therefore has a finite set of degrees of freedom, some of the 
problems must be connected with the covariance of the theory under coordinate 
transformations. The hope is of course that we are able to solve these problems 
in the much simpler model of 2+1-gravity. The big challenge will be to formulate 
a consistent second quantised theory of 2+1 -gravity and look into the problems 
of renormalisation. This issue has not been adressed to my knowledge up to now. 
Finally I would like to mention that there is a close relationship of 2+1-gravity 



with the theory of topological defects in condensed matter physics |16]. 

In this review that is based on a lecture given at the Kazan Summer School 
1995, we discuss 2+1-gravity coupled to point particles. In section 1 we shortly 
look at the solutions found in [Q and [Q. In section 2 we give the essential 
ingredients of the polygon-approach of 't Hooft. In section 3 a short introduction 
is given to the Chern-Simons formulation of 2+1-gravity. Finally we will discuss 
reference []T2[ where a new gauge is introduced which is very convenient for the 



description of particles. Almost nowhere I will go into the issues of quantisation 
deeply, but will restrict myself to general remarks. On the different ways to 



quantise the theory excellent reviews exist |[T^. I am well aware that this survey 
is only a poor selection out of a vast amount of papers that have appeared on the 
subject. I only restrict myself to the issues adressed at the summer school with 
the exception of the second part of section 1 where gravity with a cosmological 
constant is treated. 

1 The Geometrical Appraoch 

In this section we will give some simple exact solutions of the gravitational field 
surrounding point-particles. Both the case A = and A 7^ will be treated. 

In 2+1 dimensions the Einstein equations with a cosmological constant are 
given by the same expression as in 3+1 dimensions: 

G"" - Ag"'' = SnGT"" (1) 

A special feature of 2+1 dimensions is however that we can express the curvature- 
tensor in the following way: 

R''" al3 = ^'^''^^aPaGl (2) 

Using (|l]) this can be written as: 

i?'^^ „^ = SvrGe^'^V^^.rr + A{6^J^p - 5^(5^) (3) 
This implies that outside sources the curvature must be constant. 

R^R^\, = 6A (4) 
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From this we deduce that the gravitational field has no local degrees of freedom 
(there are no gravitational waves). Setting A = for the moment and considering 
a pointlike particle at the origin we have: 

= m6\r) = = z, J = 1, 2 (5) 

By symmetry arguments we have Qoi = 0. The nontrivial Einstein equations are 
given by: 

G° = -^^^^R = 87rGm6\r) (6) 

G^, = -^{D,D,-^,,D')N = (7) 

Here jij is the intrinsic 2 dimensional metric and Di is the covariant derivative 
defined with respect to that metric, ^'^^R is the 2 dimensional curvature scalar. 
Futhermore N is the well known lapse function which equals in this static case 
goo = —N'^. Taking the trace of equation (|^ we notice that we have to solve 

D'^N = DiDjN = (8) 

which is solved by N =constant. We redefine t so that = 1. In 2 dimensions 
we can always choose coordinates so that jij = e^6ij. Using this in equation (|^) 
we find: 

= -16nGmS^{f) (9) 
This is easily solved using V^lnr = 27r5^(r): 

e<^ = Cr-^^"" (10) 

Absorbing the constant C into r we finally have: 

ds^ = -de + r-^^'"(rfr2 + r'V) (H) 

Because we know that the curvature vanishes everywhere except at the particles 
positions we can transform to local fiat coordinates: 

p = — (3 = 1- AGm (12) 
e = (3^ (13) 

ds^ = -d^ + dp^ + p'^de'^ (14) 

Although the situation looks trivial now we have to be carefull. The new coordi- 
nate 6 ranges from to 2tt(3. So there is a deficit angle in space (see figure 1) The 
most important lesson we have to draw from this is that it is always possible to 
transform to coordinates in which g^y = 1]^^, but that there is a price to be paid. 
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Figure 1: The wedge cut out of space-time 



The price is that we have multivalued coordinates (or coordinates with strange 
boundary conditions). This solution is easily generahzed to N static particles. In 
this case it proves more convenient to go to complex coordinates: 

z = X + iy z = X — iy (15) 

In these coordinates the hne element for N static particles becomes: 

JV 

^ _^^2 ^ -j-j-^^ _ _ aiY^G'^idzdz (16) 

1=1 

This is a multiconical space. Transforming to flat coordinates by using the trans- 
formation: 



u= dz' n(^' - ai)"^^""' u = c.c. (17) 
•' i=i 



we construct a space from which we have to remove wedges emanating from every 
particle. In this static case it is still unimportant in which direction we choose 
the wedge. If however the particle moves, it is handy to choose the wedge behind 
the particle or in front of the particle. We will now argue why. In the static case 
the wedge is characterised by an identification rule: 

ii" = R\u^ (18) 

/ 1 

cos27ra sin27rQ; | a = AGm (19) 
— sin 2na cos 2na 



Here u"' is a point on one side of the wedge and u"" the point on the opposite side 
of the wedge. If the particle moves, the identification-rule is still the same in its 
restframe. In the moving frame we therefore have: 

= {BRB-^)\u^ (20) 



center of mass 
particle 



Hdentify; 



particle 1 



particle 2 



Figure 2: Wedge cut out around 2 particles. 

Here -B"^ is a boost matrix with arbitrary rapidity [r]) and in an arbitrary direc- 
tion. In order to avoid a time jump we have to choose the wedge symmetrically 
behind or before the particle. The effect is then that the wedge becomes a bit 
larger. This is to be expected as the energy of the moving particle is also larger. 

The particle could also carry spin (classically). The energy momentum tensor 
then looks like: 

= mS\f) = ^Se'^djS\f) = (21) 

Solving the Einstein equations is considerably more difficult {goi 7^ 0) and we will 
not repeat the derivation here. It can be found in p. Instead we will immediately 
give the result: 

ds^ = -{dt^ + 2Adtd<^ + A^d<^^) + {dr^ + r'^dip^) (22) 

Where A = AGS is the spin of the particle. The important issue is again that we 
can transform to the Minkowski line-element by the transformation: 

T = t + Aip (23) 

This has the strange consequence that we need a cut in space over which time 
jumps by an amount 8nGS. This also implies that close to the particle, closed 
time-like curves are possible. Generalisations to N massive, spinning particles 



exist |]T5[| but we will not treat that futher here. The lesson is clear: Space-times 
with N moving, massive and spinning particles can be constructed by cutting 
out wedges in space and define identifications over the these wedges. Generaly, 
these identifications are an element of the Poincare group. As an example we 
treat the 2 particle case, with total angular momentum J. We can always choose 
the particles on the x-axis. If we cut out the wedges as in figure 2 we have the 
following identification: 



u 



ai + BiRiB{\a2 - ai + B2R2B^\u - as)) (24) 
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where we didn't show the indices, and is the position of particle i. Note 
that we didn't choose the wedge of the second particle at its tail, introducing 
a time jump over the (total) wedge. The transformation is precisely a Poincare 
transformation: 

u = Au + q (25) 
A = BiRiB^^B2R2B2^ (26) 
q = BiRiB^^{a2-ai) + ai-BiRiB^^B2R2B2^a2 (27) 

If we write: 

A = -Bcom-Rcom-Bcom (28) 

this transformation really desribes a center of mass particle that is boosted to a 
speed V = tanhr^. We still have the freedom to choose the overall Lorentzframe 
so that we may take: 

A = i?com (29) 

Futhermore we may place the c.o.m.-particle at the origin by demanding = 0. 
Comparison with a spinning particle at rest in the origin suggests that the total 
angular momentum is given by the time component of the translation vector g": 

J = [5com?]° = -[BiRiB^\a2 - ai)]° = -[B2R2B^\a2 - a,)f (30) 

In the limit G ^ we can indeed recover the special relativistic result. In the next 
section we will also treat multiparticle solutions, but then we will consequently 
put the wedges behind the particles in order to make a Cauchy formulation pos- 
sible (no time jumps). The total angular momentum expresses itself then as a 
space-like translation over the total wedge. 

Next we consider space-times with A 7^ 0. The relevant reference is Again 
we consider static configurations {goi = 0) and choose conformal coordinates on a 
time=constant slice: = e'^Sij. The Einstein equations for a static configuration 
of particles (without spin) is: 



A ^ 

dd(P+-e'^ = -AnGN'^^mi5'^{z-a^) (31) 

2 i=i 

ddN + -e'^N = (32) 

dV = dV = (33) 

V = ^e-'f'dN V = ^e-'^9iV (34) 

A A 

Here we used the notation: d = 4- and d = 4=- An additional condition is 

az oz 

that initial static particles only remain static if Fqq vanishes at the location of the 
source. One can show that this implies: NV = NV = at the source. Of course 
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A>0 




Figure 3: de Sitter space-time with a particle at the origin 

we must also insist that and e*^ are real and single-valued functions. Taking 
into account the above considerations we can write down the following solutions 
for A > 0: 



N = (36) 

where f = f{z) is a holomorphic function and / = f{z) is an anti-holomorphic 
function. Also for A < we can check that we have the solutions: 



2 



|A|(1-//) 
1+ ff 



In order to reproduce delta-functions in equation (plD we demand that at z = 
we have the following singular behaviour for 0: 

~ —AGrriiN'^^ai) \yi[zz) + regular terms (39) 

The case of 1 particle sitting in the origin is now easily solved by: 

f = z-'^ f = z~'^ P = l- AGm (40) 

We can check that A^ is regular at the origin so that we can actually scale it 
to 1 at r = 0. We see that in both cases the metric becomes pathological near 
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A<0 




Figure 4: Anti de Sitter space-time with a particle at the origin 

r = 1. In the case A < this means that physical infinity is situated at coordinate 
distance r = 1. However in the case A > this is just an artifact of the coordinate 
system and we may change to coordinates that cover the whole sphere. In order 
to visualise these solutions it is convenient to use the following embedding: 

A > (41) 

U = c.c. (42) 

A < (43) 

U = c.c. (44) 

In the case A > this is just the stereographic projection, in the case A < 
this is a projection from a hyperboloid to the plane. In the /-plane the solution 
is pictured by cutting out a wedge emanating from the particle's location / = 
(as in section 1 the angle variable of / ranges from to 27r(l — AGm)). Mapping 
this on the sphere and the hyperboloid using the above embedding results in the 
above pictures (figures 3,4). 

In the case A > we see that the solution realy represents 2 antipodal par- 
ticles. In the second case we have found the Poincare-disk (Lobachevsky-space) 
from which a wedge is removed between 2 geodesies. [| It is important to notice 

^Lobachevsky, who worked and lived most of his life in Kazan, was the founder of non- 
Euchdean geometry. 
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seam 



particle 



Figure 5: Space divided up in polygons. 

that in both loop traversed around a particle would result in a transfor- 

mation: 

/ - e'^'^f (45) 

(46) 

One can check that the expressions for e"^ and are invariant under these trans- 
formations (as they should be by the requirement of single- valuedness) . In the 
case of de Sitter space one expects that static configurations of more than 2 par- 
ticles do exist. It is however not known to us if there exist explicit solutions in 
the literature. 



2 The Polygon Approach 

The relevant references are 0. The fact that the gravitational field has no degrees 
of freedom calls for an approach in which only a finite set of degrees of freedom 
survive (the only degrees of freedom of the theory are in fact the positions and 
momenta of the particles). This is precisely what happens in what we call the 
Polygon approach to 2+1-gravity coupled to point-particles. It was invented by 
't Hooft to prove that no closed time-like curves can occur in a closed universe. 
He proved that the universe would crunch before the CTC could be finished. The 
basic idea behind this method is to divide space up into rectilinear polygons with 
vertices where 3 seams meet (see figure 5) 

Just like in Regge-calculus, inside each polygon the space is flat and the metric 
is simply Minkowskian. When we move over a seam and enter a new polygon the 
coordinates we will change according to a Lorentz-transformation. Thus on every 
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seam 




Figure 6: Polygons in the x-t axis 




Figure 7: Moving seams 



polygon we choose different Lorentzframes. Futhermore we demand the following 
3 conditions: 

• On every polygon there is a rest frame such that it is an equal time surface. 

• At any time the total surface is a Cauchy surface and it is chosen to be 
equal time everywhere. 

• The metric must be continuous at the seams. 

Also particles are incorporated in this model. They sit at the end of a 1-vertex 
(see figure 5). Using the above rules and our knowledge of 1 particle solutions we 
can deduce the following rules: 

1. The lengths of the seams are equal as considered from 2 adjacent polygons. 

2. The speed of the seams is always orthogonal to the seam and equal in 
magnitude as considered from the 2 adjacent rest frames (see figure 7). 
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3. The deficit angle is always in front or behind a particle (in order to avoid 
time jumps; see figure 8). 

4. At the 3- vertices (where no mass is present), there are relations between 
the boost parameters rji and angles a, (see figure 9). These can be deduced 
from the fact that the 3 dimensional curvature must vanish at the vertex 
(although the 2 dimensional curvature may be different from zero). Roughly 
speaking; 3 quantities determine 3 other quantities. One has to take care 
however of inequalities such as the triangle inequality: |?7i| + |?72| > \rj3\. 

5. As the system starts to evolve all kinds of transitions will take place (actu- 
ally 9 different kinds). We will only picture 2 examples in figure 10. The 
transition from one situation to the next is a completely determininistic 
proces, i.e. all the new variables {rji, ai, Li) can be calculated from the old 
ones. In this way we have a completely deterministic model, free of patholo- 
gies with a finite set of degrees of freedom that can be used to study 2-1-1 
dimensional gravity. 

It is also possible to write down a Hamiltonian formulation for this model. As the 

total energy in 2+1 gravity is equal to the total deficit angle it is a logical step to 
take the total sum of all deficits as the Hamiltonian of the theory. Firstly we have 
a deficit angle of Hp = 2ap for a particle. Secondly there is a possible deficit at 
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Figure 10: 2 Possible transitions 



a 3- vertex of Hy — 27: — ay — ay — ay see figure 9) . The total Hamiltonian is 
now given by: 

H,ot^EHp, + EHv, (47) 

i j 

As configuration variables we use the lengths of the seams: Li{> 0). We already 
mentioned that is was possible to express the angles ay. in terms of 3 neighbouring 
boostparameters r)i. Also the ap^ can be expressed in terms of the mass and 
boostparameter rji at his tail. Moreover it is known how fast the lengths Lj grow, 
i.e. can also be expressed in terms of neighbouring r]i (it is constant in time). 
We can use this information together with the equation of motion: 

j^Li^{Htot,Li} (48) 

to derive that 

Pi = 2r]i (49) 
The Li, 2r]j are canonically conjugate phase space variables: 

{Li,27]j} = 5ij (50) 

The other equation of motion is trivial because the Hamiltonian H{ai{rjj)) only 
depends on the momenta: 

f^P^ = {Htot,P^}^0 (51) 

The simplicity is however a bit deceiving because we have to incorporate some 
constraints in the model. The constraints could be expected since the number of 
canonical variables exceeds the number of true degrees of freedom, which are only 
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connected with particles. Q One can very easily understand these constraints. 
First of all the angles inside the polygon must add up to (m — 2)tt where m is 
the number of angles. 

m 

Ci = ^ - (m - 2)7r = (52) 

i=l 

Secondly, the variables Lj, written as vectors, must add up to zero: 

m 

C2 = J2L^ = (53) 

i=l 

This can be reformulated in complex notation as: 

m k—1 

C2 = Y.Lke''''=0 4 = ^(vr-«,) (54) 

k=i j=i 

Futhermore we like to remind the reader that there are still the following inequal- 
ities: 

Li>0 (55) 
\Vi \ + ^ \Vk\ and cyclic permutations (56) 

which are hard to handle. The constraints (|52| , |5^ are first class and generate 
gauge transformations in the model. The first constraint evolves one polygon 
in time. As a result the surrounding seams change their position. The second 
constaint changes to a different Lorentz-frame inside the polygon. Obviously, 
again all seams start to change. It can be calculated that these changes are 
properly given by: 

6L, = {L,,Cj} 6p, = {p,,C,} J = 1,2 (57) 

One can also check that the constraint algebra closes properly. This is not a 
trivial task at all and the algebra is highly nonlinear. With that we mean that 
it is not of the simple form {Ci,Cj} = A^jCk, but on the righthand side for 
instance a sin(Cj) can appear. Of course the nontrivial part of the evolution are 
the transitions that can take place and they have to be dealt with separately. Now 
that we have formulated a Hamiltonian theory, the way to quantisation seems to 
be open. One simply changes Poisson-Brackets by commutators: 

{.,.}^-^h[.,.] (58) 

Next one chooses for instance to work in the Li representation and writes down 
a wave function: 

^d{Li,...,Ln) (59) 

•^The number of degrees of freedom for a closed (g=0) system is calculated to be AN — 12 
where N is the number of particles. 
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where D denotes its dependence on the diagram. First of all one needs to take 
into account the transitions. They must be incorporated as boundary conditions 
on the wave function. Needless to say that for a complicated diagram this is 
a very difficult task. Also the constraints and inequalities should be taken into 
account as constraints on the wavefunction to select the physical Hilbertspace. 
Although some improvement can be gained by choosing clever coordinates it still 
remains a very difficult task to complete this quantisation scheme. This will be 
one of the motivations in section 4 to go to a single-valued coordinate system. 
One can already deduce one important result in the quantum heory. As we can 
only write down expressions for sin H and cos H we only have an expression for 
the evolution-operator e~*^* and not for H itself. This implies that we cannot 
distinguish between e~*^* and Q-^iH+'iTrn)t ^ order to avoid multi-valuedness we 
have to conclude that t is integer-valued. So the model has discretised time! 



3 Gauge Theory of Gravity 

In this section we will briefly review some of the aspects of of the Chern-Simons 
approach to gravity and the way it is coupled to point-particles. The references 
are here (0)- The main result is that the Einstein-Hilbert action is equivalent 
to a Chern-Simons action where the gauge field takes values in the Poincare 
algebra IS0(2,1). If we write J'^ = |£"^'Jfec for the S0(2,l) generators and P'^ 
for the generators of translations, the algebra is: 

The gauge fields are then decomposed into: 

A^ = e/P, + ^/J, (61) 

Here " (the dreibein or tetrad) is considered as the gauge field for translations 
and 00^°" = ^e"''^'^uj^hc (the spin connection) as the gaugefield for Lorentztransfor- 
mations. The usual Chern-Simons action is: 

Ics = J d^x e^'PTriA.d^Ap + \a^A,A,) (62) 

is invariant under the following gauge transformations: 

bA^ = -V^A{x) (63) 
= -d^A-[A^,A] (64) 

The gauge parameter A can also be decomposed in terms of independent gauge 
parameters for Lorentztransformations and translations: 

A = p%x)Pa + ^%x)Ja (65) 
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The gauge transformations for the dreibein and the spinconnection are also de- 
composed into local Lorentztransformations: 

^le/ = e\,e^'K'^ (66) 

51.;/ = + (67) 

and local translations: 

52e/ = d,p'^ + e\,u;^'p'^ (68) 

52cu^" = (69) 

Witten has shown that these gauge transformations are equivalent (if the 
equations of motion are satisfied) to the usual coordinate transformations. Sub- 
stituting the expression (|6ll) into the C.S.-action and performing a 2-1-1 split of 
space-time we can write down the following action: 

/cs = j dtj ^''e. '^^^ai - Vabe^ "i^i M - VahU^o "^IN (70) 
One immediately reads from this action that s'^^e^ " is the canonically conjugate 

K \x),e^^e,\y)] = r^^'5l5\x-y) (71) 



variable to Uaf 



Futhermore, as there are no time derivatives of Cq " and uOq they act as Lagrange 
multipliers, imposing the constraints: 

Ff[cj] = e'^R}^ (curvature) (72) 

F!^[e] = e'^^ (torsion) (73) 



where: 



Rt^ = d.u^'^-d,ur + e\,u,'u^^' (74) 
T;^. = 9.e/ - 9,e, + £%,(.;, ''e/ + e/o;/) (75) 



These constraints are first class and it can be shown that they generate the 
ISO (2,1) gauge transformations: 

5G{e,uj) = {G',p(a;),Ff [a;] + Ka{x)F^\e]} (76) 

Also the constraints obey the ISO (2,1) algebra. Finally, the usual expression for 
the field-strength in a gauge theory is in this case translated to: 

F,. = [V,, VJ\ = T;^P, + RIJ, (= 0) (77) 
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The next step is to couple the particles to this action. To keep the discussion 
transparent we will couple spinless particles only, in the way proposed by Grignani 
and Nardelli 0: 



(78) 



where D^J^ is defined as the IS0{2, 1) invariant derivative: 



(79) 



This matter action is invariant under the gauge transformations 



a 



a 




(80) 
(81) 



S2q' 

52P' 



a 



a 



P 





a 



(82) 
(83) 



In the paper of Grignani and Nardelli it is also stressed that the Poincare coor- 
dinate at this stage cannot be identified with a space-time coordinate. Also 
the gaugefield " is not the dreibein or soldering form and the Poincare torsion 
is not the space-time torsion. To make a connection with a space-time interpre- 
tation we have to fix a gauge. For instance = would bring us back to the 
usual coupling of point-particles to gravity where " can be interpreted as the 
dreibein. Another possibility would however be to choose g'* = 

The C.S.-theory explained above is a theory invariant under both diffeomor- 
phisms and Poincare-gauge transformations (although they are not independent). 
One way of quantising the theory is to construct a complete set of gauge-invariant 
observables and use these as phase-space variables. In our case we must find func- 
tional of the gauge fields that are Poincare-gauge invariant and diffeomorphism 
invariant. These observables will then corresond to Hermitian operators in the 
quantised theory. In the C.S.-theory these observables can be found relatively 
easily. They are given by the Wilson-loops: 



Here R denotes the representation used for the Poincare generators. V denotes 
path ordering, 7 is a spacelike loop and Tr denotes that we have to take the trace. 
The argument of is denoted as [7] . This means that it is independent of the 
precise path of the loop but only depends on the first homotopy class of loops 
on the punctured plane. This is a difficult way of saying that all loops that can 
be deformed into each other without moving over a puncture are considered as 
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Figure 11: A Wilson- loop enclosing a particle 



the same loop. The fact that Wilson loops only depend on the homotopy class of 
loops is equivalent to the statement that Wilson loops are invariant under diffeo- 
morphisms (that are continuously connected to the identity). These coordinate 
transformations deform the path of the loop and the position of the particle in a 
continuous and invertible way: 1^^,(7) W^r(/(7))- A diffeomorphism will not 
move the path over a puncture so that it remains in the same homotopy class. 
We will now argue why the Wilson loop is invariant under these diffeomorphisms. 
It can be seen that the difference of the deformed loop and the original loop is 
again a closed loop, not containing any particles inside. Because the field strength 
inside this loop vanishes everywhere it follows from the non-abelian Stokes theo- 
rem that this Wilson-loop is actually the identity. Take for instance the simplified 
situation of figure 11: 

= WR{^i)WR{a)WR{^2) (85) 
WR{a) = TrPexpil j (fx F^[e]Pa + F^[uj]Ja] = I (86) 

It is also straightforward to show its invariance under gauge transformations 
(^?(^)): 

WR{-f) = Tr{g{xo)Vexp[fdx^A^]g{xo)-') (87) 

= TrV exp[(j^ dx>' A^] (88) 

Here we used the fact that the trace allows for cyclic permutation of matrices. 



Martin |jT9[ calculated the explicit algebra of these Wilson- variables and proposed 



to use this algebra as a starting point for quantisation. Of course in the loop 
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representation of Smolin and Rovelli it is precisely these loops that act as 
fundamental variables in the theory. The wave functionals then only depend on 
the homotopy class of loops. They proposed a transformation from the gauge 
field representation to the loop representation with precisely the Wilson-loop as 
a kernel: 

vI/([7]) = / d^i{A) Wn{\^])^{A) (89) 

where C/G is the space of all gauge inequivalent fields and djji^A) is a measure 
on this space. More about the loop representation can be found in Kirill Krasnov's 
lecture in this same volume. 



4 Gravity in 2+1 Dimensions as a Riemann- 
Hilbert Problem 

In this final section we will follow the lines of [|12| where we define a new gauge that 
proves convenient in attacking the problem of solving the gravitational field with 
point particle sources. We will work again in the A.D.M. formalism and consider 
an open universe. The hope is that we will be able to remove all redundant 
gravitational degrees of freedom from phase space by solving them in terms of 
the particles positions and momenta. This process is called reduction and works 
as follows: The total Lagrangian is schematically written as: 

^ j d'x {T^'^da,, - N^,m) + j dt Y.p,adtqt - Hm + / rf'x v^^'^i? 

(90) 

The first term is the Einstein-Hilbert action, the second is the particle action 
and the third term is the surface term needed in the case of an open universe 



1^. Our gauge choice will remove the kinetic term in E-H action. Solving the 
constraint equations at a time=constant slice will remove all Hamiltonian terms. 
Inserting the solution of these constraint equations into the boundary term will 
then generate the effective Hamiltonian. After the reduction process we end up 
with: 

/\ r r (2) 
dtpiadtql + Y^;^ j dt j d^x ^-g{g^u{p,q,x)} R{g^^{p, q, x)} (91) 

The surface term surface term is (after integration) an explicit function of p and 
q. The fact that the gravitational field carries no degrees of freedom makes it 
possible (in principle) to do this without losing any information. The problem 
of N point particles was treated in section 2 with the use of flat coordinates (the 
polygon-approach). We also mentioned that the unusual boundary conditions on 
these coordinates made them multivalued. Consider for instance a particle sitting 
at rest in the origin (see figure 12). An observer traversing via path /3i to the point 
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y-axis 




Figure 12: Particle sitting at rest in the origin 



B"^ would give this point different coordinates than an observer traversing path 
/?2- In the following we will choose coordinates in such a way that we remove the 
boundary conditions (i.e. we will have coordinates with the conventional ranges) 
but as a price we will generate a non trivial interaction term.0 If we denote by 
the flat multi- valued coordinates and by the curved, single-valued coordinates, 
the metric in terms if the becomes: 



^ u''{x)^u\x)riah 



dxf^ 



dx^ 



(92) 



First we like to choose a slicing condition. For that we view u'^{x) as embedding 
coordinates. So for every time t we want to define a function u"(x, y; t) that tells 
us how the t =constant surface is embedded in fiat 3-dimensional Minkowski 
space-time. The condition will be that locally the area of the surface must be 
maximal. This is of course always possible which assures us that the gauge will 
be accessible. Actually it is precisely the Polyakov action used in string theory 
that should be maximalised: 



Sp = j d^x ^Y'diu''dju\ 



(93) 



On the Cauchy surface we will choose conformal, complex coordinates, i.e.: 

7ij = e'^Sij z = X + iy z = x — iy (94) 
The Polyakov action then reduces to: 



d_ 

dz 



Sp = I (fz du^'du^b d=^ d = i^ (95) 



d_ 

dz 



■^This proces is very much the same as in anyon-physics where also a muhi- valued and 
single-valued gauge exist 
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It is well known that the equations that follow from this action are: 



du^du^ab = (96) 
du^du^b = (97) 
ddu" = (98) 

So u"" must be a harmonic function: u"" = u"'{z; t) + u"'{z; t). We will try to solve 
the dreibein field " = d^u"" in the following. The conditions (0,^,^) translate 
into the following conditions for the dreibeins: 

Co " = dtu" = C^{t) + dt J dz e," + dt J dz e, " (99) 

" = = "(z) (holomorphic) (100) 

" = du'' = eg "(z) (anti-holomorphic) (101) 

e.^e^^b = (null-vector) (102) 

ej^e^Vfe = (null-vector) (103) 

where C°'{t) is some vector only depending on t. Note that there is still conformal 
freedom: z f{z). These conditions in turn can be translated into conditions 
on TT*-' and jij (for the definition of the canonical momentum (see 0): 

lij = eHij TT = Tri = (104) 

The advantage of these conditions is now evident. First of all, conditions ( |104| ) 
indeed imply that the kinetic term in the Einstein-Hilbert action {n'^^ dtjij) van- 
ishes. To see that, one must split tt^^ and jij into a traceless part and the trace. 
As vr is conjugate to Tr7jj = 2e'^ we find the result. Futhermore from a math- 
ematical point of view it is very convenient that " is a holomorphic vector 
because it allows us to use the machinery of complex calculus. 

Next we will try to reduce the problem of solving the gravitational field around 
the moving particles to a mathematical problem, known as the Riemann-Hilbert 
problem. First we mention that all the information for the gravitational fields 
is really in ° and the asymptotic behaviour of the gravitational field at infin- 
ity. If we know " and the boundary conditions at infinity we can in principle 
calculate " and Cg The asymptotic behaviour can be studied by solving the 
gravitational field around a massive spinning particle in our gauge. We expect 
that at z ^ oo this is the asymptotic form of the multiparticle solution. The 
total mass of the universe is the the mass of the particle and the total angular 



momentum is its spin (see fl^). We have already argued that the flat coordinates 
are multivalued if we traverse a loop around a particle. Specifically if we move 
around one particle the result will only depend on the first homotopy class of the 
punctured plane and not on the precise path chosen. The transformation will be: 

u'^^^iBABr^-ur + qt (105) 
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where Bi is a boost-matrix, Ri a rotation-matrix and qf is the translation- 
vector, [ji] Denotes a loop in a certain class. We will look at the transformation 
properties oi e^"' — du"". 

^ {B,R,B;^ ■ e,r (106) 

The Riemann-Hilbert problem is now formulated as follows: 

Given these set of monodromy transformations (Mj = BiRiB~^) on a Rie- 
mann surface, find the vector-functions " that transform in this way if we 
move around the puncture. 

Because of the multivaluedness of " we need cuts in the plane from the 
i*'* puncture to infinity, " Must now transform as it moves over the cut. The 
familiar 1-dimensional example is of course f{z) = 2;" where / ^'^-maj 
cross the cut. To treat the R-H problem it is more convenient to use the 2 
dimensional SU(1,1) representation which is the covering of S0(2,l). The R-H 
problem still remains the same only the monodromy matrices M, now live in 
SU(1,1). Prom the spinor solution to the R-H problem (denoted by we can 
then easily reconstruct the vector solution It is also important to mention 
that the metric is single- valued (as it should be). 

g^u = Vabe,, "e, ^ ^ r)ab{Mi ■ e^Y{Mi ■ e,f = g^, (107) 

Here we used the orthogonality property of SO (2,1) matrices: 

M\M\r}ab^Vcd (108) 

Another important remark is that the monodromy around the ;2 = 00 is deter- 
mined by the other monodromies around z — af. 

= Mjv Ml (109) 

This is due to the fact that a loop around infinity (on the Riemann-sphere) is 
equivalent to a loop around all particles in the opposite direction. 

There are 2 ways of studying the R-H problem and we will only briefly sketch 
the methods here. In one approach one tries to write down a d'th order differential 
equation with N-l-1 singularities (one of which is located at infinity), d Is the 
dimension used for the representation (d=2 for SU(1,1) and d=3 for S0(2,l)) and 
N is the number of particles involved. The singularities must be of the special 
type called "regular singularities". It is always possible by a global Lorentz 
transformation to go to a frame in which the i*'^ particle is not moving. This 
then implies that the monodromy around that particle is given by a rotation. In 
this frame the spinor C° has the following singular behaviour near the particle: 

C ' {z - a,f^h^ + c^z - a,) + ...] (110) 
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where Af are called the local exponents. They must obey the following rule: 

EEA^ + EA^ = ^^^(^-i) (111) 

a=l k=l a=l 

In our case we will choose: 

Ar = (-Y,y) a, = 4Gm, (112) 

for the singularities in the finite part of the plane. The exponents at infinity are 
then determined by relations ( |109| , p.ll| ). It is important to notice that changing 



the local exponents by integers will not change the monodromy. To solve the 
R-H problem uniquely we have to choose the monodromy and the integer part of 
the exponents independently. One argument to choose certain exponents could 
be that we demand that the singular behaviour matches the behaviour found in 
pertubative calculations (small Gm). From ( [1101 ) we can easily deduce the the 
monodromy matrix: 

Mi=\^ Q j a, = AGm, (113) 

The difficulty is of course that not all monodromies commute and can be brought 
simultaneously to this diagonal form. For 2 particles however we can prove that 
the most general second order linear differential equation with 2 singularities 
and 1 at infinity can be written as a hypergeometric equation. It is long known 
that the 2 independent solutions of the hypergeometric equation transform into 
a linear combination of them if we traverse a path around a singularity. By 
matching this monodromy with the desired monodromy we can find a specific 
solution to the R-H problem 0, |5|, |l2l. 



For more than 2 particles it is convenient to to consider a (equivalent) linear, 
first order matrix differential equation of the form: 

N A. 

5r = ^F— ^ (114) 

i=i z — ai 

Y is the fundamental matrix of solutions to the differential equation. The matri- 
ces Ai are not depending on z (but may depend on Oj !). We demand that the 
solutions to this equation must fuUfiU the following conditions: 

i) Y{z,) = I 

ii) Y[z) is holomorphic in {C — {ai...aiq, oo}) 



^Recently I received a note that Grignani and NardcUi actually found this result first in the 
second reference of j?] . 
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iii) Y{z) has the following short distance behaviour near a singular point: 



Y(z) 



(115) 



where Mj = e'^'^^^' and Mj are the monodromy matrices. It is also important that 
the matrix Yi{z) is holomorphic and invertable ai z = ai. So: 



Y,{z) (r°(a,) + Y\a,){z - a,) + ....) 



(116) 



with det |y°(aj)| 7^ 0. The R-H problem is now converted into finding the dif- 
ferential equation ( |1 14| ) or equivalently to find the matrices A^. This problem 
is in principle solved in the mathematical literature |2^, Lappo-Danilevsky 



found an explicit solution to this problem in terms of a series expansion in the 
Lj-matrices. He could prove convergence of this series for small enough Lj. We 
will however not go into this technical details of this solution. Miwa, Sato and 
Jimbo found a representation of the Y{z,ai) in terms of a correlation function 
of conformal Dirac-spinors in a free field conformal field theory. The particles 
are represented by twistoperators and ensure the right monodromy properties 
if one moves around a particle. We refer the interested reader to the literature 



2e, 22, 23, 12 



5 Discussion 

In this survey we discussed some ways of treating 2+1-gravity coupled to point- 
particles. We think this is an important issue because it is the starting point for 
a quantised theory of gravity. One could for instance be interested in defining a 
consistent S-matrix for quantised particles. First of all in- and out-states should 
be defined with great care because the interaction is long range. We expect 
however that consistent in- and outgoing states can be defined because exact 
stationary scaling solutions exist (i.e. all particles move for instance towards the 
center with a constant velocity proportional to the distance from the center). It is 
however not clear if for some incoming configurations the universe crunches and no 
outgoing states exist (or there may for instance exist bound states of particles). 
Ultimately one would like to be able to describe creation and annihilation of 
particles. Concluding we would like to say that a lot of interesting work lies 
ahead. 
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